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application of these IR finite Green's functions in the case of higher order terms of the perturbative 
expansion is commented. 
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1. Introduction 



Many years ago Mandelstam obtained a solution of the SDE for the gluon propagator, in the 
Landau gauge and in the case of pure gauge QCD, that behaved as l/k 4 in the infrared [Eh. This 
solution was named as a "confined gluon solution" because it naturally leads to a linear confining 
potential. This result motivated intense phenomenological studies where the IR enhanced propaga- 
tor was substituted by more tractable functions, however these functions were still peaked at origin 
of momenta. Examples of this program can be found in the review of Ref.[Q]. 

Although not widely diffused, it is mentioned in Mandelstam's paper about the possible ex- 
istence of a massive solution for the gluon propagator, but this solution was discarded since the 
beginning in his approach. Within the same approximation ^ but with one improved three gluon 
vertex function [Q] it can be shown that this solution indeed exists in the Landau gauge. Of course, 
such soft IR behavior certainly requires a more subtle confinement explanation []5p. 

The problem with a large part of the more recent SDE calculations for the gluon propagator 
(for a review, see [^]) and the ones discussed above is that they do not lead to transverse solu- 
tions as required by gauge invariance, and it seems that the only way to have the Schwinger's 
mechanism (gauge invariant) realization in non-Abelian gauge theories is rearranging the SDE 
diagrams through the pinch technique [0, W\. This was, for the first time, devised by Cornwall 
in 1982 [Bt], who obtained a gauge invariant solution for the gluon propagator that behaves as 
\/[k 2 + m 2 (k 2 )]. As k 2 — > the function m (k ) was interpreted as a dynamical gluon mass with 
the limit m 2 (k 2 — > 0) = m 2 . In this picture the coupling constant is also IR finite with a fixed point 
behavior described by 

M0) 55 4*61n[(Lg)/A2] ' (U) 

where A = Aqcd is the QCD scale where the perturbative coupling becomes singular. Recent lattice 
QCD simulations present clear evidence for the dynamical generation of a gluon mass (a long list 
of references about these simulations can be found in Ref . [Ml and a recent simulation can be seen 



in Ref.[g]), with the SDE result for a massive gluon fitting nicely the lattice data fllOQ. 

There are several reasons to review the strong interaction phenomenological calculations on 
the light of IR finite gluon propagator and coupling constant. A simple one is that in this scheme 
we get rid of the Landau singularity in the coupling constant, which introduce singularities in 
the physical amplitudes that do not correspond to the expected physical behavior. There are also 



singularities in QCD amplitudes, like the one in the two-gluon QCD Pomeron model fllifl, that 
disappear when dealing with a dynamically generated gluon mass. The dynamical gluon mass 
also introduces a natural IR cutoff which may substitute the one that is always present in many 
perturbative QCD calculations. Moreover it is usually argued that the perturbative QCD series can 
be reorganized in order to ameliorate its behavior, but, most important of all, this optimization of 
the perturbative expansion is quite dependent on the infrared behavior of the coupling constant 



[12], and the dynamical mass generation scheme is precisely giving us one hint of which direction 
to go in order to improve the perturbative series. Finally, besides all these reasons, the fact that IR 
finite Green's functions with a dynamically generated gluon mass scale provide a better agreement 
between experiment and theory, as will be shown here, corroborates the SDE and lattice results. 
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As we shall discuss in the next section, when presenting the expressions for the gluon propa- 
gator and coupling constant, the SDE solutions cannot determine the m g value, and the best we can 
do is to determine the ratio m g /A. We will show that several perturbative calculations can be im- 
proved with the knowledge of propagator and coupling in the full range of momenta, and in some 
cases the experimental data can only be fitted with the help of these IR finite quantities. In Section 
2 we discuss that the relevant scheme to introduce these non-perturbative information into the QCD 
perturbative expansion is the one named Dynamical Perturbation Theory (DPT), then we present 
examples of phenomenological calculations that make use of IR finite Green's functions. Different 
observables are computed as a function of m g and all data indicate a small range for the dynamical 
gluon mass, what is impressive if we consider that our examples involve different hadronic mass 
scales or wave functions. In Section 3 we comment how this procedure can be extended to compute 
higher order terms of the perturbative expansion and we make a summary in Section 4. 



2. DPT at tree level 

2.1 DPT and IR finite SDE solutions 

A prescription of how the non-perturbative SDE solutions can be inserted into the perturbative 
QCD expansion was proposed by Pagels and Stokar many years ago, in the approach denominated 



DPT [13]. In their scheme the amplitudes that do not vanish to all orders in perturbation theory 



are given by their free field values, while amplitudes that vanish as X °< e ' 8 are retained, and 



possibly dealt with in an expansion in g n X. The work of Ref.[ 13] was particularly concerned with 
the effect of a dynamically generated quark mass, but as we now know from the SDE solutions that 
the gluon and coupling constant also have an infrared finite value, we can extend their formulation 
and generalize the perturbative expansion using the quark and gluon propagators with a dynamical 
mass and the IR finite coupling constant. This means that we should perform perturbation theory 
with the dressed quark and gluon propagators and the effective charge (dependent on the gluon 
mass). 



The SDE solutions and the lattice results were discussed at length in this workshop [15], there- 
fore we will not enter into details about the solutions and will just present the gluon propagator and 
coupling constant in the case that QCD generates a dynamical gluon mass, otherwise we mention 
in the references where a different SDE solution is used. We consider a gluon propagator that will 
have the form 

iA^ v (q) = P^A(q) ; P„ v = -g^v + ^ , (2.1) 

where A(q) is the gauge invariant scalar part of the gluon propagator, which in Euclidean space has 
the form 

A{Q2) ~wrkm' <22) 

The gluonic SDE solutions allow us to write a new propagator A~ 1 (Q 2 ) which absorbs all the 
renormalization group logs, exactly as happens in QED with the photon self-energy, and form the 
product d(Q 2 ) = g 2 A(Q 2 ) which is a renormalization group invariant. The dynamical gluon mass 



(mz(Q )) is given by a complicated expression falling with the momentum as 1/2 Q14p , and it can 
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be simply approximated by 



4 

171 

s Q 1 + m 1 



where m g rs ^(1.2 — 2) A, with A = Agco ~ 300 MeV. As this is still a complicated expression to 
take into account when calculating loops, in practical calculations, the best we can do is to assume 
m 2 g {Q 2 ) ?a m 2 . A simple fit for the coupling constant that is factored out in this procedure is given 

by[| 

- 2 1 

asd{q ) = 4Kb\n[{4m 2 -q 2 -ie)/A 2 \ ' (2A) 

where b = (33 — 2iif)/4Sn 2 . Eq.([2T4|) clearly shows the existence of the IR fixed-point shown in 
Eq.(l.l). It must be stressed that the fixed point shown in Eq.(|L~l|) does not depend on a specific 



process, it is uniquely obtained as we fix A and, in principle, it should be exactly determined if we 
knew how to solve QCD. The quantity d(Q 2 ) is the one that appear in all loop calculations. As 
in QED, where the only ultraviolet divergences are associated with the vacuum polarization, and 
affect the renormalization of the coupling constant, in QCD, as the pinch technique is applied and 
the vacuum polarization is summed, the only difference is that instead of a renormalized coupling 
and a massless propagator it is a factor like d{Q 2 ) that appears in the calculation, indicating the 
existence of a massive propagator and an IR finite charge. 

There are several examples of the use of DPT with the propagator and coupling constant 
discussed above [ |i"6| ] apart the ones that will be commented in the next subsections [jn|, [l8|, ^Q ], 
where the nice description of the experimental data within this procedure can be perceived. 

2.2 Pion form factor 

The asymptotic pion form factor is predicted in perturbative QCD, according to Brodsky and 



Lepage Q21|], as 



Fn{Q 2 ) = I dxj dy<j)*(y,Qv)T H (x,y,Q 2 )^(x,Q x ) , (2.5) 

JO JO 

where Q x = min(x, 1 — x)Q and Q is the 4-momentum in Euclidean space transferred by the photon. 
The function <j>(x,Q x ) is the pion wave function, that gives the amplitude for finding the quark or 
antiquark within the pion carrying the fractional momentum x or 1 — x, respectively. The func- 
tion TH(x,y,Q 2 ), is the hard-scattering amplitude that is obtained by computing the quark-photon 
scattering diagrams shown in Fig. ([I]). The perturbative kernel Th is giving by 

rp i n s 647T (2 «, t [(l-x)(l-y)«2 2 ] , l «,(xy<2 2 ) l n - 

while in the case of DPT the gluon propagator and coupling constant are exchanged by the dressed 
(non-perturbative) functions: 



T H (x,y,Q 2 ) 



~a sd (K 2 )A(K 2 ) + ^a sd (P 2 )A(P 2 ) 



(2.7) 



3 

where K 2 = (1 — x)(l — y)Q 2 and P 2 = xyQ 2 . The details of the calculation, as the dependence 



on the pion wave functions and other quantities, can be found in Ref.[17]. The comparison of the 
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Figure 1: The leading-order diagrams that contribute to the pion form factor. 



theoretical result with the experimental data can be seen in Fig.(|2j). The important point to notice 
is that this result comes out from a convolution of the gluon propagator and coupling constant with 
the pion wave functions, and, as happens in this case, for each observable that we shall comment 
in this work there will be different wave functions or scales involved in the calculation, however all 
the experimental data will indicate the same range of values for the dynamical gluon mass, which 
points to a strong phenomenological constraint on the SDE solution that we are using. Notice that 




Figure 2: Comparison of the experimental data for /^(solid line) with the theoretical determination of the 
form factor via DPT for a range of dynamical gluon masses. We also show the results obtained with the SDE 
solution of Ref.[^]. 



the purely perturbative result, obtained with the kernel giving by Eq.(2j5), does not show such a 
nice fit for the experimental data. 

2.3 Hadronic cross section in a QCD-inspired model 

According to a QCD-inspired model to compute hadronic scattering p2[], the cross section for 
producing jets with pj > pj min (through the dominant process gg — » gg) is proportional to 

o jet (s)= f dp 2 T ^f [ dx l dx 2 g(x l ,Q 2 )g(x 2 ,Q 2 ) (2.8) 

J pi dp T JX\X2>Ap\js 

min 
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where g(x, Q 2 ) is the gluon flux, and a minimum transversal momentum (pj .) establish the re- 
gion where perturbation theory can be applied. The standard calculation assumes the following 
elementary partonic cross section for gluon-gluon scattering 



d gg( s ) 



"IK 



d{s-ml) , 



(2.9) 



where mo and cCq are fitted parameters. 

DPT allow us to compute elementary cross section, like & gg (s) with IR finite quantities, where 
mo, do and p \ . are substituted by m g \ We fitted the pp and pp scattering data keeping m g as a 
free parameter. Taking a 5% variation on the minimal x 2 /DOF value indicate m g « 400^jqq MeV 
(for details, see Ref.[|rj]). One of the fits to the experimental data can be seen in Fig.(||). Again, 



even considering the quite different calculations of Eq.(2.5) and Eq.(2.8), it is amazing how the 
observables that we discussed up to now lead to the same gluon mass range, supporting the DPT 
scheme based on IR finite SDE solutions depending on a dynamical gluon mass. Moreover, the 



two different fitting parameters of Eq.(2.9) were eliminated from the calculation in our procedure! 
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Figure 3: Total cross section for pp (solid curve) and pp (dashed curve) scattering. 



2.4 Non-leptonic annihilation B meson decays 

The use of DPT can be exemplified by the diagrams of Fig.(|j), which show the different contri- 
butions for the two-body non-leptonic annihilation B meson decays in the factorization approach. 
DPT predicts that to each gluon exchange depicted in the different pictures will correspond to a 
dressed gluon propagator that enters into the amplitude calculation, as well as its coupling constant, 
i.e. the perturbative a s and gluon propagator were substituted by the product d(Q 2 ) = g 2 A(Q 2 ). 
In the standard perturbative QCD calculation we are faced with end-point divergences due to soft 
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gluon emission and an arbitrary cutoff is introduced in the amplitude: 



^=ln^(l+p^), 0<p<l. 
x A h 



(2.10) 



This divergence is eliminated in the DPT scheme with IR finite Green's functions. 




Figure 4: Perturbative diagrams leading to two-body non-leptonic annihilation B decays 



Some branching ratios of these two-body non-leptonic annihilation B meson decays are shown 
in Table ([]]). The results of Table ([[]) contain only the leading twist contribution, and we may 
expect an increase in the branching ratio values when the higher twist contributions are added. A 
discussion on the meson wave functions dependence as well as on the choice of B meson scales 
(jj, = nib or jli = nib/ 2) can be found in Ref.Jl9|]. The important point is that the best description 
of the data is giving by an IR gluon propagator and coupling constant that can be associated to a 
dynamically generated gluon mass. A full picture of these decays calculated in a totally perturbative 
approach is not so compelling as the one presented here, besides the fact that, although small, there 



is a dependence on the ad hoc cutoff discussed in Eq.(2.10). 



Table 1: Branching ratios for non-leptonic annihilation B decays obtained with the infrared finite gluon 
propagator and coupling constant discussed in the beginning of Section 2. These values were obtained with 
m g = 500MeV. A complete list of the results and experimental data can be found in Ref . [[l9|l . 



Decay channels 


m g = 500MeV 


Experiment 


&r(B° s - 


-f K+7C-) X 10 7 


1.58 


< 13.6 


£gr(B° d - 


-> K + K~) x 10 8 


7.18 


4±15±8 


@r{B° s - 


^D-k+) x 10 6 


1.54 




mr{B* - 


^D+7i-) x 10 7 


1.88 




£gr(B° d - 


^DJK + ) x 10 5 


1.98 


2.9 ±0.4 ±0.2 


®r{B- - 


^D+K-) x 10 8 


0.99 


< 1.1 x 10 5 
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2.5 A QCD-Pomeron model 



Hadronic cross sections can also be computed in the scope of Regge theory, where the ampli- 
tude for elastic proton-proton scattering at high energy is dominated by the Pomeron. Landshoff 
and Nachtmann (LN) proposed a QCD-Pomeron model where the basic Pomeron structure is rep- 
resented by two non-perturbative gluon exchange, where the gluon has a finite correlation length 



[|1 lp - The differential proton-proton cross section can be writing as 
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Figure 5: 7J and 7^. Two-gluons exchange model for the Pomeron. 



do _ \A(s,t)\ 2 

dt 16^ 2 ' K ' 



where the amplitude A(s,t) is giving by the diagrams shown in Fig.( 



A(s,t) = isSaf [Ti - T 2 ] . (2.12) 



Several results obtained for the LN model in the DPT scheme can be found in Ref.[2C] as well 



as in some of the references [16]. Actually the LN proposal match exactly with the DPT ideas, de- 
manding the introduction of an IR finite non-perturbative gluon propagator and coupling constant. 
We just show the result for the differential proton-proton elastic cross section at y/[s) = 53GeV 
depicted in Fig.@. In Fig.^ we compare different SDE solutions and the data is better explained 
by a gluon propagator and coupling constant associated to a gluon mass scale of ^(400) MeV. No- 
tice that the two gluon exchange model can explain the elastic scattering only at small t. For large 
transferred momentum it is necessary to add a three gluon exchange with opposite parity. Again 
we do have different wave functions and scales that have to be introduced in the calculation (see 



Ref.|2C]), but we still have a good agreement with the data for the same range of dynamical gluon 



masses. 



3. DPT at higher orders of the perturbative expansion 

The question that we would like to discuss in this section is that if we can have precise phe- 
nomenological tests of the a s and gluon propagator infrared behavior at higher orders of the per- 
turbative expansion. This type of test together with the ones discussed in the previous section can 
definitively provide a strong test for the SDE solutions. These SDE solutions can be applied in the 
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Figure 6: Differential pp elastic cross section at y/(s) = 53GeV computed within the Landshoff- 
Nachtmann model for the Pomeron, using different infrared couplings and gluon propagators obtained from 
DSE solutions. 



DPT scheme at the loop level, however this is not a trivial matter and requires the use of the pinch 
technique ^ to disentangle the different contributions that come from different Green's functions 
inside loops. What we will present here is a preliminary result of the DPT application in the case 
of the Bjorken sum rule and many aspects of this calculation are still in progress [j23|]. 



It has been pointed out that the Bjorken sum rule [ 24 ] could be used to know the a s behavior 



up to low energies. This is, for instance, the point of view followed in Ref.[25] using the idea of an 



effective coupling, and also of Ref.[26] in the case of analytic perturbation theory. Therefore this 
will be a perfect arena to test the behavior of the infrared quantities that we discussed up to now. 
The polarized Bjorken sum rule can be written as [ ^4| ] 



r p -"(Q 2 ) = S Bj = jf 1 dx [g p (x,Q 2 )-g1(x,Q 2 



(3.1) 



where g p (g") is the first spin structure function for the proton (neutron), which were measured 
recently at quite low Q 2 [pT]]. 

The QCD correction to this sum rule up to the fourth order in the strong coupling constant a s 
for massless particles and effective number of flavors rif = 3 is given by [ j28| ] 



r p r n (Q z ) = s Bj (Q l ) 



gA_ 

gv 
-20.21 



a s (Q 2 ) 

K 

2^ 3 



3.58 



130.0 



% 

QCs(Q 2 

K 



(3.2) 



where gA and gv are constants appearing in the nucleon beta decay. In principle the experimental 
data on the structure functions can be used to determine CC S (Q 2 ). 

In our preliminary calculation we have not considered the effect of dynamical masses inside 
the loops (we expect that at leading order this effect is small). We also neglect power corrections 
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to the Bjorken sum rule, because many calculations have shown that their effect is negligible, but 
we also believe that the power corrections will be softened by the infrared finite gluon propagator 
behavior [23]. Eq.( |3.2| ) was solved up to order a 4 in order to obtain the value of the gluon mass 



scale using the experimental data of T p ^ n {Q 2 ) obtained at Jlab [27]. We show our result in Fig.(^), 
which is fitted by Eq.(^4|) with a dynamical gluon mass value equal to = ^(300 — 400) Me V, 
that is compatible with previous phenomenological determinations of m g [|T6|]. The match between 



Jlab 2008 
Jlab 2004 
PDG 

Perturbative coupling 

■ m =380 MeV A=300 MeV 

g 

■ m =335 MeV A=250 MeV 




q [GeV] 



Figure 7: Effective coupling constant extracted from the experimental data of T'l "(Q 2 )- 

the effective charge a s d with the experimental data is impressive. We stress that any other different 
behavior for the infrared coupling constant visibly do not show such a nice fit for the experimental 
data. There are many points that still need analysis in this procedure: a) It would be interesting 
to have a formal demonstration of the DPT scheme realization in the case of the Bjorken sum 
rule through the use of the pinch technique, b) Dynamical gluon and quark masses, even if their 
effect are small inside loops, should be considered in future work. They should improve the match 
in Fig.(^) of the experimental data with the non-perturbative coupling in the 1 GeV region and 
c) Higher order corrections should be fully calculated in the DPT scenario in order to confirm 
our order of magnitude estimates. These calculations certainly will be quite difficult but they are 
necessary once we consider the good description of the experimental data shown in Fig.(^). 



4. Summary 



SDE solutions as well as lattice simulation of pure glue QCD are indicating that the gluon 



propagator and coupling constant are infrared finite. In Ref.[]10|] we can see a nice agreement 
between these non-perturbative methods, where the lattice data for the gluon propagator is fitted 
by a SDE solution associated with a dynamically generated gluon mass. We show that several 
strong interaction observables computed within perturbative QCD improved by the knowledge of 
the Green's functions, at the full range of momenta in the DPT scheme, also provide strong support 
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for the SDE and lattice results. The full scenario works so well that we can barely neglect the 
possible existence of the dynamical mass generation mechanism in QCD. 

All the examples of hadronic phenomenology that we have discussed are well described by a 
gluon propagator and strong coupling constant that are dependent on a dynamically generated gluon 
mass. No matter we deal with improved perturbative QCD calculations or QCD inspired models 
we verify that the experimental data is fitted with a dynamical gluon mass scale m g «0(2Agcz))- 
In many cases this mass scale helped us to reduce the number of arbitrary parameters in the calcu- 
lations. The tests that we performed are non-trivial in the sense that they result from the calculation 
of physical quantities where the gluon propagator or product of propagators are integrated weighted 
by different functions (involving different mass scales), and all quantities show agreement with the 
experimental data for gluon masses that are in the same range of masses predicted by Cornwall 
several years ago It is hard to believe that such coincidence is a fortuitous one. 

A preliminary account of a phenomenological test of SDE solutions at the loop level was dis- 
cussed in Section 3. The simple analysis introduced there is giving a signal that the DPT procedure, 
with the inclusion of IR finite Green's functions, is quite promising, and details about this approach 



shall be presented elsewhere Q23|]. 
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